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1. Introduction 

Since the work of Kontsevich [Q] it is well appreciated that matrix model techniques are 
very effective in solving topological string theory questions. Recently Dijkgraaf and Vafa 
proposed a new matrix model description for the topological string on non-compact Calabi- 
Yau threefolds, which are ALE fibrations of A-D-E type over the complex plane ||]. Various 
checks have been made on the claim that the matrix model solves the closed topological 
string on these backgrounds [||, || and gives the exact holomorphic information of Af=l 
quiver gauge theories, namely the superpotential and the gauge kinetic terms ||. Orbifold 
and orientifold constructions widen the range of geometries and the field theories questions 
that can be addressed, see e.g. [||, f|, ^, ||] and references therein. 
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The loop equations encode the symmetries of the matrix model as Ward identities, 
which are in principal strong enough to determine the correlators. Moreover the geometry 
of the local Calabi-Yau manifold is directly described by the large N limit of the loop 
equations for the resolvent, which encodes the spectral density of the matrix. 

One of the most elegant methods to obtain these equations is to rewrite the matrix 
model correlators in terms of CFT correlators in a theory of free chiral bosons. Correla- 
tors of products of bosonic currents are identified with matrix model correlators involving 
polynomials of the resolvent. Particular combinations of normal ordered products of the 
bosonic currents can be systematically identified with so-called Casimir fields. Their cor- 
relators fulfill a vanishing condition, which translates in the desired loop equations. This 
method is well known [10, |ll|, ^] and it was used in the context of the matrix model 
for the topological B-model to establish those terms in the loop equations of the ^4-quiver 
matrix models, which lead to the identification with the A singularities in the fibre of the 
local Calabi-Yau space Q. The Casimir fields form a W-algebra, which is the symmetry 
algebra of the matrix model [13, 14]. As opposed to a normal Lie algebra, the commutator 
of two modes in a W-algebra does not necessarily give only a finite sum of modes, but can 
also result in an infinite sum of products of modes. The special W-algebra formed by the 
Casimir fields has been called a Casimir algebra [15, 16 1. 

In the present paper we want to formulate the approach precisely enough to derive 
explicitly the exact finite N loop equations for A r and D r quiver theories, which are our 
main result and are given in equations ( |3.14| ) and ( |3.43| ) , ( |3.44| ) . The techniques of section || 
equally apply to the i?-series, but we do not work out loop equations explicitly for this case. 
The large N loop equations enable us to identify the A r and D r fibre geometries together 
with the subleading terms, which encode the resolution of the singularity by renormalisable 
deformations, i.e. complex structure moduli of the topological string B-model. 

These data calculate the exact 7V=1 gauge theory information. The finite N loop 
equations encode the exact terms of the gauge theory coupled to gravity. In the topological 
string context they encode the g > genus correlators, which was checked for Ai using the 
formalism of [17] in || for genus 1. The equivalence between the generalised Konishi 
anomalies of the M=\ theory and the loop equations, observed in the large iV-limit in [18], 
was generalised to the gravitational sector in [ 19 , [2(]] . 

The occurrence of the chiral CFT in the derivation of the loop equation seems auxiliary 
at the first glance. However it is known that the Virasoro constraints, which naturally 
appear in the CFT formulation, characterise topological string correlation functions for 
topological gravity, see [21] for review. The higher spin fields used to obtain the higher 
loop equations should lead to similar constraints for the topological string theory correlators 
in the quiver backgrounds. Recent results [22] for the open topological string established 
a chiral boson on a Riemann surface, which encodes the essential part of the Calabi-Yau 
geometry, as actually describing the complex structure deformation of the topological B- 
model. The Ward identities for the chiral boson correlators are strong enough to solve the 
all genus open string amplitudes for the topological vertex. Open/closed string duality 
extends this logic to the closed topological string. One should find a CFT formulation, 
which is identified with Kodaira-Spencer gravity for the closed string, and use the W- 
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constraints to directly fix the all genus partition function of the topological string on the 
quiver geometries. 

Another motivation to study quiver matrix models comes from |23|] . There it is argued 
that r-functions should be thought of as "the next generation of special functions" , whose 
properties should be investigated on their own right and which will be relevant to describe 



partition functions of string theory. This program was started in [23] with the r- function 
given by the hermitian one matrix model, and among the obvious next candidates would 
be the quiver matrix models. 

Acknowledgements. We would like to thank V. Kazakov, I. Kostov and C. Schweigert 
for helpful discussions. I.R. thanks the SPhT Saclay for hospitality, where part of this 
work was completed. S.C. is supported by DFG Graduiertenkolleg 271/3-02 and I.R. is 
supported by the DFG project KL1070/2-1. 

2. A-D-E quiver matrix models and free bosons 
2.1 Loop equations and matrix models 

Let us summarise the method we will use in the derivation of the quiver loop equations 



in the well-studied case of the hermitian one-matrix model, see e.g. the review [32]. The 
model is defined by the matrix integral 

d$ e -^ w W = / d\ 1 ...d\ N Y[(\ m -\ n ) 2 e-^ E ^ W(Xk) (2.1) 

m,n=l 



The first integral is over all hermitian NxN matrices, while the second integral amounts 
to expressing the first one in terms of the eigenvalues of In the notation used below, 
this is the ^4i-quiver model with the real axis as integration contour. The observables of 
the matrix model are traces of powers of <I>, 

Z mm [TY($ mi )---Tr($ m ")] = (const) Jd$ Tr($ mi ) • • • Tr($ m ") ei W($) . (2.2) 

These correlators can be conveniently expressed in terms of a generating function 

Z mm [u(zi) ■ -u(z„,)] , (2.3) 
where we introduced the resolvent uj(z) of the matrix model, which is defined as 

w( z ) = -i-Tr— *— = - V^- . (2.4) 

fc=i 

Expanding ( |2.3| ) in powers of Zf. gives as coefficients the correlators (^^). The resolvent 

obeys the so-called loop equation, obtained by a change of variables <I> — > + e/(z— 3>) in 



the integral (|2.1|) , 

Z mm [cj(z) 2 - JL-W'(z)u;(z) +p(z)] = , (2.5) 
Ng s 
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where p(z) is a polynomial of maximal degree equal to twice the degree of W'(z). The loop 
equation determines the expectation value of the resolvent in the large N limit, and it can 
also be used to set up an iterative procedure to find its l/7V-expansion [24, fT^ . 



There is an alternative way to derive the loop equation (^5|) by relating the matrix 
model to the conformal field theory of one free boson in two dimensions via Z c f t = Z mm 
JToj , |Tl| , |l2| |. Here Z c ^ is a shorthand for the vacuum correlator of the free boson in a 
complicated background which encodes the matrix integral. This relation is further used 
to map CFT correlators of n-fold normal ordered products of the free boson current J(z) 
to matrix model correlators of polynomials in the resolvent, 

Z c{t [:J(z)---J(z):] = (const') Z mm [(tj(z) - ^-W'(z)) n ] . (2.6) 

To recover the loop equations one needs in addition the relation 

dz 1 



, Z cit [T(z)] = , (2.7) 

1(2 2m z—x 

where T(z) = ^:J(z)J(z): is the holomorphic stress tensor of the free boson and C is an 
integration contour specified in section |2.6| . Also in this section, it is explained how ( |2.C| ) 
and ( |2.7D can be derived from the explicit form of Z c {t. Prom ( |2.7D one finds that Z c ft [r(z)] 
is regular on all of C, and in fact one can show that it is equal to a polynomial p(z) of 
degree less or equal to twice the degree of W'(z). Combining with ( |2.6| ) one finds 



p{z) = ±Z cft [:J(z)J(z):] = (const") Z mm [ (oj(z) - -L-W'{z)f] , (2.8) 



which is the same as the loop equation ( |2.5| ) after expressing p{z) in terms of p(z). 

For the quiver matrix models, this procedure is generalised in several ways. The quiver 
models are multi-matrix models (section ^^) and accordingly we have not only one, but 
several resolvents uji(z). Further, these resolvents obey not only a loop equation of second 
order in the uJi(z), but there is also a set of higher order loop equations. Attempting to find 
these equations by change of variables in the matrix integrals turns out to be extremely 
tedious, and exploiting the relation to CFT proves a far more efficient method. The CFT 
now consists of several free bosons (section |2.3| ), and not only the stress tensor leads to a 
loop equation, but there are also currents of spin higher than two, which obey a relation 



analogous to (2/7) and give corresponding higher order loop equations (sections |2.5| and 

U). 

These loop equations are not only essential for solving the matrix model, but describe 
also the background geometry of the type IIB string, whose low energy limit is the J\f=l 



ADE quiver gauge theory M, see section |3.3| and 3.5 for the identification of the corre 



sponding geometries from the large N loop equations. 
2.2 The A-D-E quiver matrix model 

In general, a quiver matrix integral of rank r takes the form [25, g] 

Z mm = (const) y II d< 2™ e -^ W( *' Q) . (2.9) 

i=l (m,n) 
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Here the $j are Ni x Ni matrices and the Q mn are N m x N n matrices. The potential 
W($>, Q) is given by 

r r 

w($, Q) = s v Qn ®j Qoi + Yl w ^ ' ( 2 - 10 ) 

i,j=l i=l 

for some polynomials Wi{x). The constants s^,- are antisymmetric, Sj,- = — sji, they obey 
Sjj = 1 if i < j and the nodes in the quiver diagram are linked, and = otherwise. 
The notation (m, n) for the range of the product in ( |2.9| ) denotes all pairs (m, n) with 
1 < m,n < r s.t. s mn 7^ 0. The definition of the integration region requires some care, see 
also (27]. We will address this problem in section 2.4. 

In principle one can proceed by integrating out the Q mn and expressing Z mm as an 
integral over the eigenvalues Xij of ^j, where / = l,...,iVj. This would result in the 
expression 

„ r Ni r N, Nj 

z mm = / i[d\ k>K j] n (\u - A^) 2 n n n - a jv)~ m e- s , (2.n) 

J k,K i=l I,J=l i,j=l 1=1 J=l 

KJ i<j 

where 

S = lj2Wi(k, L ) , (2.12) 

9s l,L 

and the constant in ( |2.9[) has to be chosen appropriately. In practice this would require a 
consistent definition of the integration regions for the ctnd. Qmn 

to avoid divergences. 

For the purpose of this paper we will take ( |2.11| ) - or rather a regularised form thereof, 



described in section 2.4 - as a definition of the matrix model and think of ( |2.9D as a 
motivation for considering an integral of the form ( |2.11| ), 

From the general quiver models one obtains the A-D-E matrix models by choosing the 
|sjj| to take the special form 

= 2(5jj Aij , (2.13) 

where Aij is the Cartan matrix of a rank r Lie algebra of A-D-E type. In this case we have 
A^ = ctj), where the on are the simple roots. The integral ( p. 11 ) can now be written 
more compactly as [25, [10|, |ll|, |26|, 

„ r Ni Nj 

Z mm = / l[dX k>K H H [J(Ai,i - Xj,j)^' a ^ e- s , (2.14) 

k,K i,j=l 1=1 J=l 

where one has to choose an ordering of the pairs (i, I). Here we have chosen (i, I) < (j, J) 
if either i < j or else if i = j and I < J. Different choices of ordering change equation 
( |2l4|) by a factor of ±1. 
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have [/(l)-charge q and conformal weight h = \q 2 ■ The OPEs involving J(z) read 

1 a 

J(z)J(w) = -tj + reg(z-w) , J{z)V q {w) = —^—VJw) + reg(z-w) (2.15) 



2.3 A free boson representation 

In this section we will construct the integral ( 2. 14] ) as a correlator in the chiral CFT of r 
free bosons JlO| , [Tl] , |l2| ]. To fix conventions, consider one free boson p(z) for a start. 

The [/(l)-current idip(z) will be denoted by J(z). The vertex operators corresponding 
to the chiral part of the exponentials :e tqifi ^ : of the free boson field are called V q (z). They 

1 

(z—w) 2 ' ov ' v ' yv ' z—w 

and the conformal block with n insertions on the complex plane, together with a charge q 
placed at infinity, is given by 

n 

WM-vMW = s qi+ ... +qn>q J] (zi-zjy^. (2.i6) 

i,i=i ; i<j 

Here (q\ is an out-state dual to the highest weight state \q) of C/(l)-charge q. By definition 
we have 

(q\V qi ( Zl )...V qn (z n )\0) = lim Li 2 (0\V- q (L)V qi ( Zl ) ■ ■ ■ V qn (z n )\0) . (2.17) 

L^oo 

The modes of the current J(z) are defined via 

J(z) = Y,Jk z-*- 1 , J k =i ^z k J(z) , (2.18) 



where 70 is a circular contour around the origin. The OPEs ( p.!5[ ) imply the commutation 
relations 

[Jm, Jn] = m5 m+nfi , [J m , V q (z)} = qz m V q (z) . (2.19) 
Using these it follows that for k > 

n 

J k V qi (zi)---V qn (zn)\0) = {Y,^) k )V qi { Zl )-..V qn {z n m , (2.20) 

i=l 

and thus also, for any out-state (X\ and k > 0, 

(X\e- tJ >°V q M)---V qn (z n )\0) = e-^M* (X\V qi ( Zl ) • • -V qn (zn)\0) . (2.21) 

Let us now return to the general case. Denote by g a simply laced Lie algebra of rank 
r. Let K be the Killing form on g and fix once and for all a basis H a , a = 1, . . . , r of the 
Cartan subalgebra h of g s.t. K(H a , H b ) = 5 a ^. Let h* be the dual of h and denote by e a 
the basis dual to H a . Another basis of h* is provided by the r simple roots Oj of g. The 
bilinear form on h* induced by K is denoted by (•,•). 

Consider the chiral CFT consisting of the product of r free bosons, with components 
j( a \z) and V q a \z), where a = l,...,r. To an element u = ^ a u a e a of h* assign the 
current J u (z) and the vertex operator V u (z) as 

r r 
r{z) = Y j n a J {a) {z) and V u {z) = \{V^{z) . (2.22) 

a=l a=l 
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One computes the conformal weight of V u (z) to be h(V u ) = \{u, u). With the conventions 
(2.22), the analogues of formulas (2.15), ( |2.16| ) and fl2.19|) are given by 



J u (z) J"(w) = + ieg(z-w) , J u {z)V q {w) = { -^V q (w) + ieg(z-w) 



{z—w) e 
(q\V qi ( Zl )---V gn (zn)\0) 



J 9iH \-q n ,q 



n 



(Zi - z. 



z—w 



(2.23) 



[J^J v n ] = {u,v)m8 m+nfl , [J^,V q (z)} = (u,q)z m V q (z) . 



Here u, v, q as well as q\, . . . , q n are elements of h*. The modes of J u (z) are written as J%. 

To the r potentials Wi(x) in ( |2.10| ) we assign a h*-valued function W(x) as follows. 
Suppose 



W i (x) = J2t 



for i = l, 



,r 



(2.24) 



k=0 



and define vectors S h* via (r^aj) = t$. We then set W(x) = Ylk>o T k xk e h*. This 
potential enters the definition of an operator H, which encodes the potentials of the matrix 
model in the free boson representation, 



9 s o=l ^70 Z7 ™ 5s fc>Q 



As a final ingredient we need the screening charges 

dz 



Qi 



2m 



V ai (z) , 



(2.25) 



(2.26) 



where are some integration contours, to be specified in section 2.4. Note that because g 
is simply laced we have (ai,cti) = 2 for all i and the vertex operators V ai (z) indeed have 
conformal weight one. 

The free boson representation of the quiver matrix integral Z mm in ( 2.14 ) is now given 



by 



11 



Z c{t = (N\e~ H e Ql ■ ■ ■ e Qr \ti) where (N\ = (N iai +- ■ -+N r a r 
which is equal to Z mra up to a iVj-dependent constant 



■'eft 



Off Z n 



■ N 



where C 



1 



N 



Ni\---N r \ 



(2.27) 



(2.28) 



To see this equality note that 



1 



J cft 



JVi! • - - JV r ! 



(N\e 



k (Qi) Nl ---(Q r ) Nr \o) 



1 



Nil ■ --Nrl 



n 



dX 



k,K 



2-ni 



e b. 



(2.29) 



k,K J ~< k "" j,J 

where in the first step we inserted the definition of H and used that by conservation of the r 
[/(l)-charges only one term of the exponentials of the screening charges in ( 2.27 ) can con- 
tribute. In the second step the screening charges have been replaced by the corresponding 
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Figure 1: The original integration contour 7 and the shifted and truncated contours 7fcV used in 
the regularisation prescription. 



integrals and the relation (2.21) has been employed to remove the modes JT k . Inserting 



the explicit expression for the r-boson conformal blocks in ( p. 23 ) one recovers the integral 



(2.14). 



2.4 Integration contours 

As already mentioned, in this paper we will consider (|2.9|) as a formal expression motivat- 



ing the eigenvalue representation (2.14), which we will take as a definition of the A-D-E 
quiver matrix model. However, even in (2.14) the integrand contains potential singular- 
ities, namely whenever (aj,ay) < and two eigenvalues belonging to the nodes i and j 



approach each other. Below, in a procedure similar to [27], the expression ( 2.14 ) is defined 
as a limit of a regularised integral. The regularisation involves a number of ad hoc choices, 
but these turn out not to affect the loop equations we will derive. 

Let 7 : M -> C be a smooth contour parametrised such that it has unit tangent 
|7(i)| = 1. Given a small e > define a family of contours jl'x on the interval [L, —L] as 



T&(*)=7(*) + 



I £ 



K + £ k r} N. 



■7(0 , 



(2.30) 



'3=1 



see figure |T[ With these contours we define a regularised integral J re s( e ' L ) [ • ] as 



I reg ( £ ' L )[/(A 



r N k 

nn 



IK=1 J ^k,K 



2m 



-)/iA,. :-... ) • 



(2.31) 



Here it is implied that the contours jt'ic are integrated on the interval [L, —L] only. 

Using the integral operator ( 2.31| ) one defines the regularised partition functions Z^^ £ ' L ^ 



and Z™± ■ With this choice of contour the various eigenvalues A/% k are always at a finite 



v cft 

distance from each other and Zmm £,i ' 



^reg(e,£)y£, - g & ^g m -(- e ^ numD er, because it is 



defined as a multiple integral over a finite region of a bounded function 
Below we will be working with 



lim Z reg(e ' L) 

r mm 



(2.32) 



which a priory still depends on e. However the loop equations turn out to be independent 
of this parameter. Note also that requiring the limit ( 2.32j ) to be finite imposes further 
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constraints on the initial contour 7. The regularised integral ( 2.14j ) converges if (but not 
only if) Re(Wi(-f(L))) -^ooasi^ ±00 for all l<i<r. 

2.5 Higher spin currents commuting with screening charges 

The loop equations in the matrix model will be constructed from fields in the chiral algebra, 
which commute with the screening charges. These fields will be called Casimir fields for a 
reason explained in section 

Before turning to the definition of Casimir fields we recall some notations to express 
the OPE of chiral fields. We will use the conventions of 28], section 6.5. For two chiral 
fields A(z), B{z) the chiral fields {AB} n (z) are defined via 



{AB} n (w) 



A(z)B(w) = £ 7" JW ;7 . (2.33) 
^— ' [z— w) n 



n=— 00 



The (generalised) normal ordered product (AB)(z) of two chiral fields is defined to be 

, . „x , s I dx A(x)B(z) , , _ , . . 
{AB){Z) = % 2ri x-z = {ABUZ) ' (2 - 34) 

where 7 Z is a circular contour winding tightly around z. Multiple normal ordered products 
are defined recursively, e.g. {ABCD){z) = (A(B(CD)))(z). In general this form of normal 
ordering is neither associative nor commutative, see e.g. fE8| appendix 6.C, or the appendix 



of 1 15] for details. We will also be using this generalised normal ordering for free boson 
fields, where it reduces to the usual notion of moving creation operators to the left. 
Consider a general spin s chiral field W(z), 

s 

W ( z ) = Yl Yl Yj d Z\:.Z n {d mi J (ai) ■■■d m "J^){z) . (2.35) 

n=l ai,...,a n mi,...,m„>0 

miH \-m n =s— n 

We will give two equivalent definitions for W(z) to be a Casimir field. 

Definition 1: The spin s chiral field W{z) is a Casimir field if VoH^(0)|0) = 0, for all 
V{z) = V ai (z), i = 1, . . . ,r. Here Vq denotes the zero mode of the current V{z). 

Definition 2: W(z) is a Casimir field if, for i = 1, . . . , r, the OPE of V ai (x) with W(z) is 
a total derivative in x, 

V ai {x)W(z) = ^-(^(x-z) n A n (z)) , (2.36) 

n 

where the A n (z) are some chiral fields. 

To see that the two definitions are equivalent, first note that, by the state- field corre- 
spondence, VbW(0)|0) = is equivalent to § j^zV( x )W(z) = {VW}i(z) = 0. Inserting 
( |2.36 ) in the last expression, we see that definition 2 implies definition 1. Conversely, if 
{V r H^}i(z) = we can write 

v(xW(z) _ v ggnW _ a / A ivwUA ) (2 . 37) 

v ' v ; ^ (x-zY dx\ ^ 1-n (x-z)™- 1 / v ; 

n=— 00 x ' n=— 00 v ' v ' 
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Figure 2: If a Casimir field W(z) is analytically continued through a screening charge integral, the 
screening integral can be deformed back to its original form. Casimir fields commute with screening 
integrals. 



so that definition 1 implies definition 2. 

As we have seen while arguing the equality of the two definitions, a Casimir field W(z) 
has the property <f ^V{x)W{z) = 0, i.e. a contour integral of V(x) around W(z) can 
be left out. It follows that W(z) does not have a discontinuity across the screening charge 
integration contours, as illustrated in figure [2|. This is the first important property of 
Casimir fields. The second property appears when integrating W(z) around a regularised 
screening charge contour 7^V as in figure ||, where 71 = jt' K and yo = 71 yi = 71 (L). 

The function f(z) is assumed to be analytic inside the z-integration contour 72, and will 
later on be set to f(z) = (z — for some point x outside the contour 72. The double 
integration indicated in figure |3| can be solved 

/ 2^ / ^n f{z) ' ^ = A{VI) ~ A{m) ' (2 - 38) 



with A(y) given by 

^»)-E E yy^fW ttW,,). (2.39) 

fc=0 «=o v ' 

In verifying this calculation one can substitute the OPE ( |2.37D , then Taylor-expand both 
{yVF} n (z) and f(z) around the point y and carry out the z- integration, which gives rise 
to a Kronecker delta, removing the n-summation. The y-integration can be solved trivially 
because of the derivative d/dy introduced by the OPE. 
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The ability ( |2.38 ) of W(^) to replace screening charge integrations by field insertions 
at the integration boundaries is crucial in the following argument, which ultimately leads 
to the loop equations. 

2.6 From Casimir fields to loop equations 

As a generalisation of ( |2,27 ) let us introduce the notation 

^cft |>i(*l) •••&»(*»*)] = (iV| e - H e^---e Q ^i(x 1 )--^ n (x n )|0} (2.40) 

for the (unnormalised) expectation value of some chiral fields <pk{z) in the matrix model 
background. Consider the equation 

&jf ssi^^'M")]- - < 2 - 41) 

u lev 

where the contour 7 ev surrounds all contours 7tV, but not the point x. This equation was 
given in [jD], 14, 10] in terms of modes, and in integral form (with the L-limit implicit) in 



e 



1 26, [j^] (for the stress tensor) as well as in ||. Equation ( |2.41| ) will be established in detail 
in appendix [A]. The outline of the argument is as follows. The integration contour 7 CV is 
deformed to encircle each one of the regularised contours 7^'^-. By ( |2.38|) such an integral 
can be replaced by an insertion of A(y{) — A{dq). In the limit L — > oo this insertion will 
cause an exponential damping in the correlator, due to the non-trivial out-state (iV|e -H . 

The rewriting of ( p. 41 ) as a loop equation will proceed in several steps. First note the 
commutation relations 

[Jk, H n ] = -k 0(—k) (u, T^ k ) n H™- 1 , (2.42) 
9s 

which can be verified by recursion. Here 6{x) is the Heaviside function. As a short calcu- 
lation shows, this in turn implies that 

H J u (z) = (j u (z) - -(u,W(z)))e- H with W(z) = ^-W(z) . (2.43) 

^ 9s ' oz 

Define the positive and negative mode part of J u (z) as 

J«(z) = ^ Jk z'^ 1 and J^(z) =J2 J k z ~ k ~ l ( 2 - 44 ) 
fc>0 fc<0 

so that J u (z) = Jl{z) + JH(z). Since (N\J% = for all k < it is easy to see that 
(N\J u (z) = (N\ J+(z). In fact, using the definition of the normal ordered product one can 
verify the relation 

{N\(d mi J Ul ■■■d m "J u ")(z) = (N\d mi Jl 1 (z)---d mn Jl n (z) . (2.45) 
Putting all this together we arrive at the conclusion 
z & [{d mi J Ul ■ ■ ■ d mi J Un ){z)] 

= (N\ ]J (d mk Jl k {z) (u k ,d mh+1 \N(z))) e- H e Ql ---e Q ''|0) 

k=l 9s 



(2.46) 
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1 1 

Wfc ( z ) = nY1 ~T\ where N = N 1 + --- + N r 



On the matrix model side, we have r different resolvents Uk(z), defined as 

•V ^ z A/, A 

K=l K ' n 

Expectation values of the resolvents are defined as 

Zmm[w fll (zi) . ..U aa (z 8 )] 



(2.47) 



jreg 



r Ni Nj 

■ ai ( Z1 )...u aa (z s ) n n \v x i - e-t^i 
i,j=i i=i j=i 



1 Ei, l Wi(\l) 



(2.48) 



(i,I)<(3,J) 

The link between ( |2.46 ) and ( 2.48 ) is provided by the equations 



E 



and 



9s > - 

fc>0 



N ^ 



-k-i d 



(a) 



dr., 



,-fc-l 9 



J^\z)e- H 



k>0 



dt 



Here we have set rjf^ = (t*., e a ). Recall that by definition t^' = (r^, a«) = X^a=i( a «> e a) T fc"' '> 



Ui(z)e~ 



(0 



(2.49) 



(2.50) 



which results in a corresponding relation for the partial derivatives, namely 

8 E' ' a 

i=i 



9r ; 



(a) 



ft; 



dt 



(2.51) 



In this way one obtains, for example, 



2™ft^[jW(z)] 



fc>0 



-i d 



(a) 



I(e a ,W'(z))) Z 
5s 



reg(e,L) 
eft 



iV^(e a ,a 4 )^^z 



i=i 



-fc-i 9 



fc>0 



(0 



C^N Z^ L) [ ^(e a ,a t )co t (z) 



l(e a ,W(z))) CffZSt^ 
9s ' 

' 'e«,W'(*)) 



(2.52) 



8=1 



Ng. 



Note, however, that in writing (|2.50| ) we have expanded [z— \h,K) = z Ylj>o(^k,K / ' z) 3 . 
This is possible only if \z\ > \Xk k\- For finite L we can always choose a large enough Rq{L) 
such that all contours 7^'# are contained in a disc of radius Rq{L) centered at the origin. 
From ( |2.49| ) and ( 2.50| ) we thus know that ( 2.52 ) holds for all z outside of this disc, i.e. 
\z\ > Rq(L). For \z\ < Rq(L) the relation then holds by analytic continuation in z. Since 
we have established ( |2.52| ) for all values of L and z, it will continue to hold if we take the 
limit L — ► oo. 

In general we set, for u G h*, 



Vu{z) = ^2{u,ai)uJi(z) 



i=l 



Ng s 



u 



W(*)) 



(2.53) 
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In the L — > oo limit we then obtain the following relation between the CFT expectation 
value of a normal ordered product of free boson currents and the matrix model expectation 
value of a polynomial in the resolvents: 

(d^J^---d m "J u -)(z) ] = C^N n Z^[d^y ul (z)-..d m -y Un (z) ] . (2.54) 

This equation, together with ( 2,35[ ) and ( p. 41 ) gives the desired relation between a spin s 
Casimir field W(2) and a loop equation of maximal order s in the resolvents. Concretely, 
suppose we are given a Casimir field in the form ( |2.35| ). Then 

f ^nzh Z ^ (£ ' L) = C " N ' {F(X) + P(X)) ' (2 - 55) 



where 



F(x) = -Y J N n - s £ d " n 

n=l ai,...,a a m\,...,m n >0 

mi-\ \-m n =s—n 

x Z™t> L) [ d^y {ai) (x) ■ ■ ■ d m -y (an) (x) ] (2.56) 



P(x) = Y J N n - s E d > 

n=l ai,...,a a mi,...,m„>0 

m\-\ \-m n =s—n 

dz 1 



mi...m« 
ai...a„ 



2ni z—x 

Too 



Z^> L )[d™y {ai) (z)---d^y M (z) ] (2.57) 



as can be seen by deforming the contour j ev into a sum of two contours, one encircling the 
point x and the other, denoted by 700, a large circle containing all the eigenvalues as well 
as the point x. Here we also abbreviated yr a )(z) = y ea (z). 

The function P{x) is analytic in the entire complex plane. Furthermore the lhs of 
(2.55) behaves as c/x for some constant c as x tends to infinity. On the other hand, since 
for the resolvents we have u>i(x) = Ni/N ■ x _1 as x — > 00 and since the potential terms 
with higher derivatives are subleading in x, the function F(x) behaves as 

F{X) = J^ Z ^' L) £ <1K. W 'W)-K. W 'W) +0(x sD ~ 1 ), (2.58) 
where 

D = max( deg(W/(x)) \ i = l,...,r) . (2.59) 

In order for F(x) + P(x) to behave as c/x for x — > 00, P(x) thus has to be a polynomial 
of maximal degree sD. Taking the L — > 00 limit in fl2.55| ), one finally arrives at the loop 
equation associated to a Casimir field in the form ( |2.35| ), 



, 



«[5> n ~ s E E C-.C"5 mi y {ai )(x)...^ yK) (x) + P s (x) 

n=l a\,...,a a rai,...,m„>0 

WilH \-m n =s—n 

(2.60) 

where P s (x) = —P{x)/Z mm is a polynomial of degree deg(P s ) < sD. 
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2.7 Undetermined parameters in the loop equations 

In order to use the loop equations to determine the matrix model correlators recursively, 
we need to know more about the a priori undetermined polynomials P s (x) appearing in 



(2.60). Let us for this section assume that all of the Wi{x) have degree D+\. Writing out 
the product of y's in the integrand of P(x) in ( |2.57| ), one finds the leading terms to be of 
the form (W'(z)) s + lo(z)(W (z)) 3 ' 1 + ATW '(z)(W (z)) s ~ 2 + ..., where we have omitted 
all sums, indices, etc. Because of the asymptotics of u>(z) one sees that the (a priori 
unknown) resolvents start to enter P(x) only at order x^^" 1 ' 1 , while the coefficients 
of x Ds to x- ^ -1 ) are directly determined through the potential W(x). It follows that the 
polynomial P s (x) entering a loop equation of order s has D(s—1) undetermined coefficients. 

On general grounds one expects an independent loop equation of order s for every Lie 
algebra Casimir of g of order s. This will be motivated in section [|, where the relation 
between Lie algebra Casimirs and a special kind of W-algebra, called Casimir algebras, is 
discussed. 

The Lie algebra g has r Casimir operators, whose order s is related to the exponents e of 
g via s = e+1. The exponents obey the relation X^eecxp(g) e = R+> wnere R+ = \ (dim g—r) 
denotes the number of positive roots of g. The overall number of undetermined coefficients 
in the polynomials P s {x) of the r independent loop equations now reads 

#(undet. coeff) = ^ D(e+1-1) = DR + . (2.61) 

eSexp(g) 

The same number of free parameters appears when setting up a large iV-expansion. 
The eigenvalues of the matrices are located on the cuts of the resolvents iVi(z). For small 
enough g s the cuts will be located around the solutions of the classical equations of motion, 



which correspond to the critical points of the quiver potential ( 2.10 ). There are DR + such 
critical points (2^, and for the large N expansion one has to fix the corresponding 

DR+ filling fractions. One can thus expect a (however quite nontrivial) bijection between 
the filling fractions and the parameters Q2.61|) . 

It should be remarked that the large N expansion is in general just a formal procedure 
and does not necessarily reproduce the value of the matrix integrals for partition function 
and correlators, see [31]. 



2.8 The quadratic loop equation as an example 

In this section we recover the familiar feature that the quadratic loop equation is linked to 
the Virasoro algebra 13, 14], see also e.g. [11, 32 1. 



Recall that the OPE of a weight h primary field <j>{z ) with the conformal stress tensor 
T{z) is of the form 

7 \ Q 

T(z)<j)(w) = rj0H + — <t>(w)+reg(z-w). (2.62) 

[z— w) z z—w ow 

This can always be written as a total derivative in w if the conformal weight h is equal to 
one. Since all the screening charge currents V ai (z) are primary and have conformal weight 
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one, T(z) is a Casimir field. In terms of the [/(l)-currents of the r free bosons we have 

T(z) = j2\(J {a) J {a) )(z) . (2.63) 



2 

0=1 



Applying equation ( |2.60| ) to this Casimir field and substituting ( 2.53 ) yields the quadratic 
loop equation for the quiver models 

Z%*[ ^AijUiWufa) - -^j^W'ix)^) + P(x)] =0 (2.64) 

i,j=l 9 s i=l 

Here Ay = (aj, otj) is the Cartan matrix of g and we used the identity (ojj, W(x)) = Wi(x), 
which follows from the definition of W(x) below ( 2.24 ). Thesummand (W'(x), \N'(x))/(Ng s ) 2 



has been combined with 2P2(x) to the polynomial P(x), which has maximal degree 2D, 
but at most D coefficients which are not directly determined by W(x), see section 



Quadratic loop equations have been found for similar multi- matrix models in 26, 12 1, and 



in |33[ the large N limit of (2.64) was given for the ^4 r -case. 
3. Examples 

3.1 A r quiver model loop equations at finite N 

In this section we want to apply the techniques introduced above to find a set of r loop 
equations for the ^4 r -quiver matrix model. We have 

exponents of A r = {1,2,..., r— 1, r} (3-1) 

and thus expect the loop equations to have orders 2, 3, . . . , r, r+1. We will start their con- 
struction by introducing a convenient set of r+1 vectors {s\, . . . , £ r +i} in h*, corresponding 
to the weights of the vector representation of A r , 

r r k 

£ i = 5^ afc ~ ^2 ak ' for i = 1, . . . ,r+l , (3.2) 

k=i k=l 

where a\, . . . , a r are the r simple roots of A r . Using also (c^, aj) = A$j = 25ij — (Jij+i — 
Sij-i one can verify the properties 



ai = Ei-e i+ i , ^£fe = , (e i ,£j) = 5i > j — . (3.3) 

k=\ r 

As a further ingredient we will need the commutation relations of the modes JV; with 



the modes of the currents V ak (z), which we will denote by V Qfejm . These follow from ( 2.2S| ) 
to be 

[J^, V ak!n ] = (u, a k )V ak>m +n , where V a ^ m = j> ^^^ki 2 ) ■ ( 3 - 4 ) 

In writing the contour integral for the modes V ah m we used that V ak {z) has conformal 
weight one. Using (|3.3|) and (|3.4|) , a short calculation shows that, for any /i€C, 



[V ak> o,(» + Jl\)(» + f k 1 +1 ))=J« l {V ak ^ 1 -V ak> - 2 . (3.5) 
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Proceeding analogously as in |34], 16] we define a state \R r +i) (closely connected to a 
quantum Miura transformation), 



r+l 



\Rr + i) = Y\(n + Jl\)\0) , (3.6) 



1=1 



which will serve as a generating function for r Casimir fields. 

In fact, consider a field W(z) of the form ( 2.35| ), but without derivatives. Writing out 
the definition of the normal ordered products, one verifies 



W(0)\0)= d ai ... aa J^ ) ---J { L{ ) \0) ■ (3.7) 

ai,...,a s 



In the same way, the field R r+ i(z) corresponding to ( |3.6[) , i.e. the field with the property 
i? r +i(0)|0) = |-Rr+i) is just given by the normal ordered product 

R r+ i(z) = (A 1 A 2 . . . A r+1 ){z) where A k {z) = /xl + J £h {z) . (3.8) 

We would like to show that, for any /jgC, R r +i(z) is a Casimir field. To this end we 
verify definition 1 given in section 2.5 . Consider the transformations 

Vct k ,o\Rr+l) 

fc-1 r+l 

=U(fi+ji\)[v ak , ,(»+ji\)(»+f\ +i )] n (m+w) 

i=l i=k+2 
fc-1 r+l 

=n^+ j -i)( j -^^-i-^.-2) n (^+ j -i)i°> (3 - 9) 

i=l i=k+2 
r+l 



n (m+ ^om^-i-^^io). 



i=l 
i^{fc,fc+l} 

Here, in the first step we used that by (|3.4|) Kt fei o has nontrivial commutation relations only 
with Jt\ and J_* +1 . In the second step (3.5) is substituted, and finally this part of the 
expression is commuted past the right product, using that V akyTn has trivial commutation 
relations with Jf^ if i > A;+2. 

Next we will investigate the states J"^V^ fc! _x|0) and V^ fc) _2|0), and show that they are 
in fact equal. This will imply that Q3.9| ) is equal to zero, and hence R r+ \(z) is indeed a 
Casimir field. 

First note that by definition 

\a k ) = V ak (0)\0)= j —-V ah (z)\0) = V akt - 1 \0) . (3.10) 

Further, the state V^ fc] _2|0) has Lo-eigenvalue 2 and charge However, all states with 
Lg-eigenvalue 2 and charge are of the form J" 1 |afc), for u € h*. To find which is the 
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correct value of u, we compose both states with (a k \Jf from the left. Using ( 2.23| ) and 
fl3.4| ) one easily checks 



(a k \JiJli\a k } = (v,u){a k \a k ) and {a k \JiV ak _ 2 |0) = (v,a k ){a k \a k ) . (3.11) 

For this to hold true for any v G h* one needs to have u = a k . We have thus established 

V ak ,- 2 \0) = r\\a k ) . (3.12) 

Knowing that R r +i(z) is a Casimir field for any \x G C, we can expand in and find the 
individual Casimir fields (z) of spin s = 1, . . . , r+1 as, see ]34j] and also e.g. |35|, [jl| [l6| , 

W ( - S \z)= E {J £ ^ ■ ■ ■ J £l °){z) . (3.13) 

l<ii<-<i s <r-+l 

Note that = since the e k sum to zero. We can now apply the relation ( 2.6Q ) 

between Casimir fields and loop equations to all of the W^ s \z). The result can again be 
conveniently written in terms of a generating function 

r+1 r+1 

2£&[ n(M+w*)) + e^ 1- *^ ] = ° . ( 3 - 14 ) 

A:=l fc=0 

where P k {x) is a polynomial of degree < kD. The functions y £k (x) are expressed in terms 
of the resolvents and the potential according to ( 2.53j ). Setting uiq{x) = = uj t+ i{x) we 
can write 

y £k (x) = Lo k (x) - u k -i(x) + t k (x) . (3.15) 

Using (|3.2| ) together with (aj,W(x)) = W!(x) we find that the tree level potentials of 
quiver model action ( |2.12| ) enter in the combination 

ys ys i=i i= k 

Note that ( ^To]) can easily be inverted, 

= E(W X )-^( X )) • (3-17) 

Expanding ( 3. 14] ) in [i gives rise to r loop equations for the A r -quiver matrix model, which 
hold at finite N. The equations arising at the powers and fi r are trivial. 

3.2 A closer look at the cubic and quartic loop equations 

The cubic loop equation of the A r model corresponds to the coefficient of /x r_2 in the 
//-expansion of ( p. 14 ). A priori it is given as an ordered sum over y £i y £j y £h . The ordered 
sum can however be simplified by rewriting it in terms of sums over the full index range 
and using YH=i Vei = 0. Explicitly 

^ r+1 ^ r+1 - r+1 

E VnVejVei. = 3| E VeiVeiV*k ~ 2 E W^e* + 3 E^^ ' ( 3 - 18 ) 

l<i<j<k<r+l i,j,k=l i,k=l i=l 
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where on the rhs all but the last term sum to zero. This yields the following form for the 
cubic loop equation, valid for general r, 



. r+1 



z reg 



[E(^) 3 + 3P 3] =0- (3.19) 



mm | 

" i=l 



Then, using u>o = = u> r+ i, Ej = + £j + i and defining r\ = Ng s , one obtains 

T 

i=i 11 11 (3.20) 

- -WiW/)(e i+ i,W')) + Qsl =0, 
77 r] /J 

for some polynomial Qs(a;) of degree < 3D. This form of the cubic can be transformed 
further by making use of the quadratic loop equation ( |2.64j ). Substituting £j + i = £2 — 
Y^\=2 a k i n to ( |3.2C| ), the coefficient of (£2, W) is just a polynomial (due to the quadratic 
loop equation) and can be absorbed into Q3, resulting in 

r ^ ^ 

^mm[E (wi-lW? -W f 2 _ lWi + -W/(u, 2 - -W!tJi 

1=1 V V i (3-21) 

~ E K 2 - ^^+1 - ^to) E w * + ^ 

' i=2 ' fc=2 

where Qs(x) is some polynomial of degree < 3D. For r = 2 this equation reduces to the 
one in J?]]. 

The quartic loop equation can be treated in the same way, taking the coefficient of /j, r ~ 3 
in ( 3.14| ). For ^3 this coefficient is just y ei ye 2 y £3 ye A and the loop equation reads explicitly 

^mm [ (3Wi + 2W' 2 + W^- Ar]U X ){W[ - 2W' 2 - W£ - 4r/cji + At]lo 2 ) 

(W{ + 2W' 2 + 3W£ - 4r/w 3 )(W / i' + 1W' 2 -W^- Ar]uj2 + 4r/w 3 ) (3.22) 
- 256r/ 4 P 4 ] = . 

For general r it is convenient to rewrite the ordered sum over four indices similarly as was 
done above for three indices, giving the A r quartic loop equation in the form 

r+1 r+1 

2S&[ EWW-2EW 4 + 8P 4]=0. (3.23) 

i,k=l i=l 

3.3 A r quiver model loop equations at large N 

In the large N limit (keeping Ni/N and g s N fixed) the ra-point functions in the matrix 
model factorise into products of n one-point functions. The loop equations encoded in 
( 3.14 ) become algebraic. Abbreviating y £k (x) = ZS[i/ £k (i)]/fi and replacing /i-tzwe 
can write 

r+1 r+1 

H(z+y £k (x)) = -E^(*K +1 - fc > (3-24) 

k=l k=0 
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where again Pk(x) is a polynomial of degree < kD. It follows that, given the polynomials 
Pk(x), the y £k {x) are minus the r+l roots of the polynomial in z on the rhs of ( p. 24 ), i.e. 
minus the r+l solution to the equation 



z r+l 



Y,Pr+i-k(x)z k = , (3.25) 



k=0 



where we also used that ( |3.24| ) forces Pq(x) = — 1 and P\{x) = 0. This also proves a claim in 
p| that the loop equations of the ^4 r -quiver matrix models reproduce, in the large N limit, 
a deformed reduction to one complex dimension of the singular, non-compact Calabi-Yau 
threefold geometry associated to the ^4 r -quiver. Explicitly, the non-compact Calabi-Yau 
three-fold is given by the fibration of an ^-singularity over the complex x-plane, [36, 29, 2|, 



r+l 



u 2 + v 2 + Yl(z + t k (x)) = . (3.26) 



k=l 



Here the tk(x) as given by ( p.!6| ) are polynomials, fixed by the tree level potentials, and 
hence the geometry reduced with respect to the u and v direction is a nodal curve. In the 
quantum geometry as described by the matrix model loop equation the tk(x) are replaced 



by y tk (x) and the DR + parameters dicussed in section 2/7 are turned on. The curve does 
not factorise algebraically any more and the double points get resolved. Note that the 
number of parameters DR + , which are not fixed at tree level is sufficient to resolve all 
double points. 

Topological string amplitudes of the B-model and exact gauge theory quantities are 
calculated in terms of periods of the resolved curve with respect to the meromorphic dif- 
ferentials, which are reductions of the Calabi-Yau (3,0) form over r cycles of ( |3.26| ). They 
are given by 

r ^ 

r]i = y ai (x)dx = [J2 A ij^j(x)--—W^x))dx, (3.27) 
j=i 9 s 

where we set u = aij in ( PD and used W'(x)) = W((x). 

Note that expanding out ( |3.24D gives the r loop equations as ordered sums, 

E ye n (x)...y £ Jx)+P s (x) = . (3.28) 

l<ii<— <i s <r+l 

A corresponding equation has been conjectured in [^] to arise from an analysis of the 
Konishi anomaly in quiver gauge theories. Using Yli=i Ve^x) = 0, equation ( j3.28| ) can be 
rewritten in the more concise form 

r+l 

i=i 



for some polynomial Q s (x) of degree < sD. The equivalence of ( 3,28| ) and ( |3.29|) can 



be seen recursively. One first checks the statement for s=2. For general s one rewrites 
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the ordered sum in terms of several sums over the full index range and uses the already 
established identities ( 3.29[ ) of degree less than s. For example in the case s=3 one rewrites 
the ordered sum as in ( 3.18| ), where only the last term survives. In the case s=4 the 
same procedure leads to the large N form of ( 3,23j ); using the case s=2, one can replace 
£3=iy £l (z)%(z) 2 = Q2(x) 2 to obtain the form ( p9j) . 

3.4 D T quiver model loop equations at finite N 

Next we turn to the investigation of .D-series. For the Lie algebra D r with r > 3 we have 

exponents of D r = { r— 1 , 1, 3, . . . , 2r— 3 } . (3.30) 

Accordingly we expect a generating set of Casimir fields of spin r and spins 2,4,..., 2r— 2. 
The procedure of |p7|| , which we will recall below, is to start by constructing the Casimir 
field W^(z) of spin r and then find the fields of spin 2,4,... 2r— 2 as elements of the OPE 

W i - r \z)W i - r \w). 

As for A r we first need to choose a convenient set of vectors £{ in h*. In the case of D r 
we choose r such vectors, which form a basis of h*. Let the simple roots of D r be numbered 
such that the Cartan matrix Au = ctj) is of the form 

Aij = 25ij — 5ij + i — 5ij-i + 5i jr (5j jr -i—5j iT -2) + 5j >r (Si >r -i— <5i, r -2) • (3.31) 

The vectors Ej, for i = 1, . . . ,r are defined as 

r-2 

Si = Y, ak + IK + (-l) 5i ' r "r-l) , (3.32) 

k=i 

where for negative range the sum is taken to be zero. The vectors ei,...,e r have the 
properties 

(£j,£j) = 5ij , Ui = Ei — (1— S ijr )£i + l + <5j, r £ r -l • (3.33) 

Using the relations ( |3.4| ) and ( |3.12| ) (which are not specific to A r but valid for any simply 
laced Lie algebra) one can check the equations, for i < r— 1, 

v ai)0 Ji\JT\o) = J-\v ai ,-i\o) - v ai ,- 2 \o) = o , 

(3.34) 

VcrfiJ-TJ-M = -J-iV ar ,-i\o) + v ar ^ 2 \o) = o . 

Define the field W( r \z) via 

W (r) (0)\0) = JV-i • • • J -i\°) ■ ( 3 - 35 ) 

Since all of the r J-modes entering (|3.35| ) commute and also each of the screening charge 
zero modes V ait o commutes with all but two of the J Sk one finds 

r 

^,0^(0)10) = ( [J Ji\)v ai , Ji\Jir\o) = o , 

k=l 

(3.36) 

r-2 

V ar , W^(0)\0) = ( H Jl^Va^Jlr'Jl^O) = . 

k=l 
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Thus (z) is indeed a Casimir field. Since by translation invariance also [V^o, (z)] = 
0, we have in particular that, for i = 1, . . . , r and all values of z, 

V aifi W^(z)W^(0)\0) = . (3.37) 

It follows that upon writing out the OPE 



2r 



W (r) (z)W {r) (0)|0) = z"" 1 T^ (2r - m) (0)|0) 



(3.38) 



every coefficient iy( 2r - m )(0)|0) is a state corresponding to a Casimir field. This gives an 
a priori infinite number of Casimir fields, but they will not all be independent. However, 
one can show that the W^ r \z) together with W^ s \z), s = 2, 4, ... , 2r— 2 do indeed form 
an independent set of Casimir fields, see [37| and e.g. |l6| section 6.3.3. 

Let us investigate the fields W^ s \z), s = 2, 4, . . . , 2r— 2 in more detail. In terms of the 
modes W m of W^(z), defined via W^(z) = J2m z ~ m ~ rw m, the OPE ( ggg ) takes the 
form 

W^(z)W^ (0)|0) = z~ m ~ r W m W- r \0) , (3.39) 



so that W^(0)|0) = W r - a W- r \Q). Using the definition of the normal ordered product, it 
is not difficult to show that 



W<r\z) = (^■■■J^)(z) = £ z-(™+-+™r+r) :J^...JX: 



W„ 



E 



mi,...,m r S 



• J mj 



r 



(3.40) 



mi,...,m r GZ 
mi-\ \-m r =m 



where : • • • : denotes the usual normal ordering of modes, such that the positive modes are 
to the right. Because of ( 2.23 ) and ( |3.33 ) all the J-modes in ( 3.40 ) commute and the 
normal ordering can be omitted. If a term in the sum ( |3.40D for W m is to contribute to the 
product W m PF_ r |0), we need the to be in the set S = {1} UZ<_i. This gives explicitly 



mi,...,m r £S 
miH \-m r =r— s 

s/2 

= E E 

n=l l<«i<---<j n <r mi,...,m„>0 

miH \-m n =s—2n 



(3.41) 



The second expression for W"^ s ^(0)|0) is obtained form the first by cancelling all Jf.* 
where = 1 and redefining the other as — m^— 1. 

To proceed we will express the higher J-modes ( |3.41| ) in terms of derivatives. Noting 
that (L_i) n J_i|0) = n! J_ n __i|0) we see that the field corresponding to the state J_ n _i|0) 
is ^d n J(z). Using this, we finally find 



s/2 



w'-H*) = E E 



E 



( J £l i d mi J £ 'i • • • J £ ^ d mn J £l « ) (z) 



n=l l<i\<--<i n <r mi,...,m n >0 

miH \-Wm=s— 2n 



777-1 ! • • • m n ! 



(3.42) 
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With the help of ( |2,60| ), we can now write down the r loop equations resulting from the 
Casimir fields {W^{z), (z), (z), . . . , W^ 2r - 2 \z)}. For (z) one finds 

Z^[y ei (x) yea {x) ■ ■ -y £r (x) + P r {x) ] = , (3.43) 

where P r (x) is a polynomial of degree < rD, while for W^ s \z) one gets 



~ ljH:, 1 i.r)<)""!h„{.r) 

m\-\ \-m n =s— 2n 



E^ E EII k m J k + ] = > ( 3 - 44 ) 

n=l l<«i<-<i n <r mi,...,m„>0 fc=l 



where P a (x) is a polynomial of degree < sD. The y ek {x) are related to the resolvents via 
( 2.53 ). Explicitly, defining also u>q(x) = 0, 

Ue k (x) = LOk{x) - Uk-\{x) + 5 k) r-lUr(x) + t k {x) , 



r-2 



t k (x) = --L(e k ,W'(x)) = --L(j2Wl(x) + \{W r {x) + (-1)^(3 
iV 5 S iV 5 S V 

Again this relation can be inverted to give the resolvents Wj in terms of the y e ., 



(3.45) 



1 * 1 

Uiix) = -(2-5 i>r -i-5i, r ) E {ve k {x) - t k {xfj - -5 itT -i{y £r {x) - t k (x)^ . (3.46) 

fe=i 



Note that in contrast to A r , the -D r -loop equations (3.44) in general contain derivatives 
of the variables y £k {x). These derivatives appear because the corresponding Casimir field 
( 3.41 ) contains modes J_ m with m > 1. Recall that this was not the case in the generating 
function (|3.6| ) for the A r Casimir fields. 

However, three of the Casimir fields can be expressed through the modes alone. 
The first is, by definition, W^ r '(z) in (3.35). The second and the third are the stress tensor 
T{z) and a spin 4 field U(z), defined as 



no)|o) = ±E J Viilo) = ^ (2) (o)|o) , 



i=l 



[/(o)|o} = ± e ^i^i^i^iio)-iE J -i J -i J -i J -ii°) 

l<£<7<r i=l 

= ^( 4 )(0)|0) -L_ 2 L_ 2 |0) . 



(3.47) 



The last equality can be verified using 

L ^EE«-* : • ( 3 - 48 ) 



2 

fcez i=i 



Further, to see that U(z) is indeed a Casimir field one can use that W^\z) is a Casimir 
field and that the screening charge zero modes commute with the Virasoro generators, 
[Kii,0)im] = 0. This in turn follows from the commutation relation of the L m with the 
modes K n of ctny primary spin one field K(z^), which read. ^L rnj -^-n] — — ^-^-"n+m" 
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The Casimir fields W^ r \z), T(z) and U(z) thus result in loop equations which do 
not involve derivatives. Note that this fits well with the observation that D3 = A3, as Lie 
algebras, and the three Casimir fields of ^3 with spins 2, 3 and 4 do not contain derivatives. 
For the W^ s \z) with s > 6, linear combinations similar as in ( |3.47 ) do not seem to exist. 



3.5 D r quiver model loop equations at large N 

When taking the large N limit of (3.44) only the coefficient of N° survives. For this 
coefficient we have n = s/2, so that it does not contain any derivatives. Abbreviating 
again y £k (x) = Zmm[y £fc ( x )]/-^mm the large TV limit of the loop equations ( |3,43| ) and ( 3,4-J| ) 



reads 



and 



Veii.x)y £2 {x) ■ ■ ■ y £r (x) 



P r {x) 



ys n {xf ■■■y eia/2 {xf 



(3.49) 
(3.50) 



Ps(x) , 

l<H<---<i s /2<r 

where s = 2, 4, . . . , 2r — 2. As opposed to their finite N form, the large N loop equations 
can easily be written in terms of a generating function, 



\\{z + y £k {xf) = -Y,P2r-2k{ 



X )Z 



(3.51) 



fe=i 



fc=0 



where we set Pq(x) = — 1 and p2 r {x) = —P r (x) 2 . Note that ( p. 51 ) is weaker than the 



set of equations (|3.49D and ( 3.50| ) because the coefficient of z° just gives the square of 



relation (|3T49| ). From (|375lD we see that the y £k (x) are equal to ±i times the r zeros of the 
polynomial in z given by X^fc=o P2r-2k{x)z k = with the constraint that the signs of the 
y £k (x) have to be chosen such that equation ( 3.4S| ) holds. 

Similar as in the A r case, the large TV" loop equations of the D r -quiver matrix model 
is the deformation of the reduction of a Calabi-Yau three-fold to a nodal curve. The non- 
compact threefold is a fibration of a -D r -singularity over the complex x-plane with defining 



equation [36, |29| , 



u 



r r r 

+ v 2 z+ -(Y\(z + t k {x) 2 ) - Yl t k (x) 2 ) +vH t k (x) = , 



(3.52) 



k=l 



k=l 



k=l 



where t k {x) are given by ( |3.45| ). The deformation of nodes is achieved by by replacing t k {x) 
with y 6k (x) and turning on the DR + parameters. Note that due to the Z2 symmetry of 



( 3.51 ) pairs of douple points are resolved by this. The reduction of the holomorphic (3,0) 
on the curve gives meromorphic forms as in ( p. 27 ). 



4. Relation to Casimir algebras 

The aim of this section is to relate the calculations in sections [2] and | to a special form 
of W-algebras, the so-called Casimir algebras. We will start by some general comments on 
W- algebras. 
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4.1 Some generalities on W-algebras 



By a chiral algebra one denotes a subsector of a full conformal field theory which consists 
only of holomorphic fields, i.e. of fields 4>(z,z) which obey d/dz <f>(z,z) = 0. The chiral 
algebra need not contain all such fields, the only requirements are that it closes under the 
OPE and that it contains the conformal stress tensor T(z). 

Chiral algebras can be studied as objects on their own right, without reference to the 
original full CFT. This leads to the notion of a conformal vertex algebra, see e.g. for a 
mathematical exposition of the subject. In fact, a standard approach to study CFTs is to 
start with a chiral algebra and to use its representation theory to construct the full CFT. 
This has lead for example to the Virasoro minimal models and the WZW-models, see e.g. 
& 

For the Virasoro minimal models, the chiral algebra is just the Virasoro algebra, gener- 
ated only by the stress tensor T(z) itself. For WZW models the chiral algebra is a current 
algebra, that is, it is generated by fields of spin one. An obvious generalisation is to study 
chiral algebras that are still generated by a finite number of fields, but where the fields can 
have any integer spin. These are the W-algebras; an extensive review can be found in pq|. 



W-algebras made their appearance in [39], where a chiral algebra generated by T(z) and 
an additional field of spin three was investigated. 

A qualitative distinction between the Virasoro or current algebras and general W- 
algebras is that the OPE of two generating fields is no longer a linear expression in the 
generators. On the level of modes this implies that the algebra of modes of the generating 
fields is not a Lie algebra. The commutator of two modes can contain (infinite sums of) 
products of modes. 

W-algebras appear in the context of matrix models and integrable hierarchies in the 
form of W-constraints. In this case one has a set of differential operators W n which anni- 
hilate the r-function of the hierarchy and which do not form a Lie-algebra. Instead their 
commutators can give non- linear combinations of the W n , and the resulting structure is 



that of the mode algebra of a conformal W-algebra |L^, 14]. 

A link between the W-constraints in the matrix model context and the W-algebras of 
conformal field theory is provided by the free boson representation of the matrix model, 
see [|H], [ll|, |12|. In this paper we worked in the free boson representation right from the 
start. Another link, this time more generally between the Hirota bilinear equations for 
r-functions and W-algebras, is given by the orbit construction, see e.g. [40] and references 
therein. 



4.2 Casimir algebras and WZW-models 



The W-algebras we are interested in for the purposes of this paper are of a special type, 



called Casimir algebras. These are defined as follows, see [15] as well as [16] and references 
therein. Consider a WZW-model at level k for a Lie algebra g. The modes of the 
spin one currents T l (z) that generate the chiral algebra span an affine Lie algebra at level 
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k, which we will also denote by g^, 

dim g 

[Ti,Ti}= f\Tl +n + kmK{T\Ti)5 m+nfi , (4.1) 

l=\ 

where T % denote the generators of the underlying finite Lie algebra g, the f 1 ^ are the struc- 
ture constants of g and K(-, •) is the Killing form on g. Note that g is canonically embedded 
into gfc via T % \— > Tq. Denote this embedding by i. Similar to the coset construction of 
conformal field theory we can consider all elements in g^ that commute with all elements 
of t(g). 

Definition: The Casimir algebra W[gfc/g] is defined as the (vertex-) subalgebra of the 
chiral algebra g^ obtained by taking all fields 4>(z) of g& that obey [Tq,</>(z)] = for all 
generators T % of g. 

Via the state- field correspondence we can give an equivalent characterisation of W [g& / g] 
in terms of states in the vacuum module of g^, 

[T*,<K*)]=0 & T V(0)|0)=0. (4.2) 

The chiral algebra g& is spanned as a vector space by fields of the form 

W^ mn (z) = (d mi r* • • • d m "T")(z) . (4.3) 
One can verify (see |]16|]) that, for generic level k, the linear combination 

W mi .., mn (z)= dii-in^XW ( 4 - 4 ) 



'1 



is in W[gfc/g] if and only if £^ i ^ii...i n ^ 11 ■ ■ ■ T l " is in the centre of the universal en- 
veloping algebra of g, i.e. if it is a Casimir element of U(g). This is the reason for the 
terminology 'Casimir algebra'. For certain specific values of k, however, one looses the 
'only if in the above relation (just take the free boson stress tensor in the free boson 
realisation of gi as an example) . 

4.3 Casimir fields and Casimir algebras 

We would like to find the relation between the notion of a 'Casimir field' used in sections 
|2| and || and the Casimir algebras introduced above. 

Choose a Cartan-Weyl basis for g. In the free field realisation of gi, the free boson 
currents J^^z) provide the gi-fields H l (z), where H l denotes an element of the Cartan 
subalgebra of g. The gi-fields E ai (z), for the ladder operators E° l E g and oti a simple 
root, are given by the free boson vertex operators V ai (z), up to a cocycle factor, which we 



do not spell out explicitly, see e.g. 28] section 15.6.3. We would like to establish 



<p(z) is a Casimir field <^ <fi(z) G W[gi/g] . (4.5) 

It is enough to check the defining condition in Q4.2; ) for the modes Eq ' and Hq for i = 
1, . . . ,r. If (j)(z) is a Casimir field, then jQ a ^(0)|0) = since by construction the state 
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0(0) 1 0) is in the vacuum module of the £/(l) r -algebra and thus has C/(l) r -charge zero. 
Furthermore, by definition 1 in section |2,5| we have 1^00(0)10) = 0, so that indeed 4>{z) S 
W[gi/g]. 

In the converse direction we have to be more careful, because in the definition we used, 
Casimir fields are build only from the gi-currents H l (z), corresponding to the Cartan 
subalgebra of g, while the fields in W[gi/g] are constructed from all generators T l {z) in 
gi. Suppose 4>(z) G W[gi/g]. Then because Hl<fr(0)\0) = we know that the state 0(O)|O> 
is indeed in the vacuum module of the C/(l) r -algebra and can thus be expressed in terms 
of the modes alone. The remaining conditions Sq'0(O)[O) = then imply that 4>{z) 
is a Casimir field. 

In the g = A r and g = D r examples discussed in section || we have constructed r fields 
of W[gi/g]. One can now wonder if by repeatedly taking normal ordered products, these 
fields (and their derivatives) generate all of W[gi/g]. There are arguments (see |16|] and 
references therein), using the Drinfeld-Sokolov reduction and the character technique, that 
this in indeed the case. 



A. Appendix: Integral form of W-constraints 

In this appendix we establish equation ( 2.41| ). Denote by I{L) the integral in ( |2.41| ) before 
taking the limit. In the correlator Z^^^[W(z)] there are N = N± + • • • + N r screening 
integrals, one along each of the contours 7^'^-. Let us denote 

e,L - f dx 



Qk(% K )=f E , L ^- V - k (x) • (A.l) 
The integration contour 7 ev of W(z) can be deformed to encircle each of the 7fc'^f- 

As 

in ( 2.38 ), carrying out the z and the corresponding screening integration results in the 
insertion of a pair of fields A(yi)—A(yo) at the endpoints of the encircled contour 7^. 
The function I{L) can thus be rewritten as follows, 

= £ ^ -h (Af|e ~ H n<w# ^)i°> 

= E^l e_H TiQAlfj) [Mv k x K ) -A k {yi K ))\0) . (A ' 2) 

k,K j,J^k,K 

Here = 7^(L), y ' = L) and A k {y) denotes the field ( [2.39| ) with V replaced 

by V ak and f(y) = (y—x)^ 1 . The individual field insertions A k (y) are a sum of terms of 
the form 

(const) • j^—^ 8<{V ak W} r (y) . (A.3) 

The state \tp) = d q {V ak W} r (0)\0} is an element in the C/(l) r -highest weight representation 
Tt ak of charge a k and with highest weight vector The space H ak is spanned by the 

states J" 1 ^.! • • • J^ n mn _i\a k ), for m, > and m G h*. Thus the field d q {V ak W} r (y) is 
itself a linear combination of fields of the form 

4(V) = (d m 'J^---d m -J u -V ak )(y) . (A.4) 
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In order to find an estimate for the integral I(L) we will first establish the formula 
(N\e-"V qi ( Xl )---V qn (xMy)\0) 



P(xi, ...,x n ,y) 



(N\e- H V qi (x 1 )---V qn (x n )V ak (y)\0) 



(A.5) 



Q(xx, ...,x n ,y) 

where Y2i 1i + a k = N, and P, Q are polynomials, whose detailed form will not matter. 
Take a field 4>{y) which is also of the form (A. 4) and consider the function g(y) defined as 

g(y) = (N\e~" V qi ( Xl ) ■ ■ ■ V qn (x n ) (d m J u 4>)(y)\0) . (A.6) 

Using the definition of the normal ordered product we can write 

g{y) = <f £p —<iv| e - H v qi ( Xl ) ■ ■ ■ v qn (x n ) d m j u {z) fa) |o) 

L^-ti,^^^^ v * M - v * i * ) j ° (z) *» m < a - 7 ' 



(N\e- H V qi (x 1 )---V qn (x n ) 0(y)|O). 



In the first step, the integration contour 7^ has been deformed to encircle the points Xi as 
well as the point 00. Also, by partial integration the m derivatives d/dz have been shifted 
from J u (z) to the factor (z— The last step uses the OPE (|2.23|) as well as the fact 
that due to ( |2.43| ) the asymptotic behaviour of the correlator in the second line of ( A. 7 ) is 
given by 



(N\e- H V qi (x 1 )---V qn (x n ) J u (z) 0(2/) |0> 



.9, 



-(u,W'(*)) (^~ H V qi ( Xl )---V qn (x n ) 4>(y)\0) + O(z- 1 ) . 



(Ai 



Furthermore, it was used that for the contour integration around infinity we have, for some 
polynomial p(z), 



dz 
2tH 



(p(z)+0(z- 1 ))(-ir 



0' 



dz m z—y 



d m p(y) . 



(A.9) 



Applying this procedure recursively allows us to strip a field (f>(y) of the form (A.4) of all 
its components d m \J Ui , establishing equation ( |A.5| ). 

With the help of ( |A.5| ) we can now find an estimate for a correlator involving screening 
integrals. Consider the equalities 

(N\e-" HQ^fj) 0(y)|o) 

d^j,J\ P(^l,U ■ ■ ■ ,K,N r ,y) /r rt, -H TTi; /l \ if / \ ln\ 
j ~2Td) Q(\,j,. ■ ■ , \ r ,N r ,y) V ^ Vm (A.10) 



n 



n 



2ni J Q(\ 1U ...,\ Nr , y ) 



e 9 
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Here P, Q, P, Q are polynomials whose arguments consist of y as well as of all A's except 
for Xk,K- The point y depends on L and is taken to have either of two values 



y = y{L) 



(A.ll) 



In the first step in ( A.1C| ) the screening integrals have been written out explicitly and 
was substituted. In the second step the explicit form of the correlator, obtained from 



(|^2ll) , ( ^23|) and ( gjg ), has been inserted and the rational part of the expression has been 
absorbed by redefining the polynomials P and Q appropriately. 

Denote by F(L) the multiple integral inside the square brackets in the last line of 
equation ( A.10[ ). For finite L this integral is finite, since the only singularities of the 
integrand could occur at coinciding integration variables Xij = Xj t j which is impossible 
by construction of the regularised contours "fi'f . Furthermore, because of the asymptotic 



behaviour of the contour 7(2:) imposed at the end of section 2.4, the integrand of F(L) 
receives an exponential damping as its arguments approach infinity. We can conclude that 
the L — > 00 limit of F(L) is well defined, 



lim F(L) = C 



L— »oo 



(A.12) 



for some constant C. For the correlator in the first line of ( A.10| ) this implies that 

. (A.13) 



hm(AV H HQihtf) <M2/)|0) = limine 



j-W k {y(L)) 



L— >oo 



Since from ( |A.2 ) we know that I(L) can be written as a linear combination of terms of the 
form ( A.lCp , the above equation proves that indeed lim^oo I(L) = 0, as claimed. 
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